Path Integrals, and 
Classical and Quantum Constraints* 

John R. Klauder 
Department of Physics and Department of Mathematics 
University of Florida, Gainesville, FL 32611 

Abstract 

Systems with constraints pose problems when they are quantized. 
Moreover, the Dirac procedure of quantization prior to reduction is 
preferred. The projection operator method of quantization, which 
can be most conveniently described by coherent state path integrals, 
enables one to directly impose a regularized form of the quantum 
constraints. This procedure also overcomes conventional difficulties 
with normalization and second class constraints that invalidate con- 
ventional Dirac constraint quantization procedures. 

1 Introduction 

In order to discuss the quantization of systems with constraints it is first 
important to briefly review what are two absolutely essential features of the 
very process of quantization itself. First of all, we hold it self evident that: 

1. The abstract operator formulation of quantum mechanics is correct 
and fundamental. 

As a corollary of this viewpoint we next observe that: 

2. In order to properly describe quantum mechanics, it is necessary that 
any functional representation of quantum mechanics have an associated un- 
derlying operator formulation. 

*Based on a presentation at the 8th International Conference on Path Integrals, PI2005, 
Prague, Czech Republic, June, 2005. 
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In particular this second property applies to: 

a) The Schrodinger partial differential equation formulation of quantum 
mechanics, 

and to: 

b) Any version of a path integral formulation of quantum mechanics. 

Although this paper is concerned with path integrals, it is pedagogi- 
cally useful to spend a few paragraphs on how these principles apply to 
the Schrodinger equation. The operator form of this equation is given (in 
units where ^ = 1) by 

id^(t)/dt^n{t)^(t) , 

where ^{t) denotes the time dependent abstract vector G i^, the abstract 
Hilbert space, and H{t) denotes the (possibly) time-dependent, self-adjoint 
Hamiltonian operator. As an example suppose that the system in question 
is a certain anharmonic oscillator characterized by the fact that 

where P and Q denote abstract, irreducible, self-adjoint Heisenberg operators 
that satisfy not only the Heisenberg commutation relation [Q, P] — il, but 
they also satisfy the Weyl form of these relations, namely that 

for all real c-numbers p and q. It was shown by von Neumann [1] that, apart 
from unitary equivalence, there is only one realization of the operators P and 

Q, namely, the Schrodinger representation P ^ —i d/dx and Q x, acting 
on the Hilbert space i^^(R) of functions ip{x), x G R. Substitution of this 
representation for P and Q into the abstract operator form for 7i yields the 
usual Schrodinger equation for this example, namely 

idij{x,f)/dt = -lip"{x,t) + |a;2^(a;,t) + AxV^^^.t) • 

So much for the obvious associations that apply to the Schrodinger equation. 

However, it is useful to enquire what may happen if the connection to the 
operator formalism is broken. The classical Hamiltonian for the anharmonic 
oscillator is normally taken to be 
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but after a canonical coordinate transformation of a suitable kind it is pos- 
sible to express the classical Hamiltonian for the same system in the form 

H{p, q) = , 
or in still other coordinates in the form 

etc. All of these functionally unequal forms properly describe the same phys- 
ical system in the indicated canonical coordinates. Clearly, to promote the 
coordinates in these distinct cases to canonical Hciscnberg operators would 
lead to Hamiltonian operators with quite different spectra and they all can 
not be physically correct. How is one to know which set of canonical coordi- 
nates to promote to canonical operators so as to obtain the correct physical 
spectrum for the specific system under consideration? The answer, accord- 
ing to Heisenberg, Schrodinger, and Dirac [2], is that the classical canonical 
coordinates should be chosen as "Cartesian coordinates" . 

To put some further substance in this remark, it is useful to appeal to 
coherent states [3]. In particular, let |0) denote a normalized vector that 
satisfies {Q + iP)\0) — 0, namely, |0) is the ground state of an harmonic 
oscillator with unit angular frequency. Let 

for all (p, q) G R^, denote the set of coherent states. Then, in view of 
the Hciscnberg commutation relation, it follows for a general Hamiltonian 
operator 7Y(P, Q) that 

H{p,q) = {p,q\n{P,Q) \p,q) 

= {0\n{P + p,Q + q)\0) 
^n{p,q) + 0{h;p,q) ; 

the last form of this expression is particularly evident for polynomial Hamil- 
tonians. In any case, in the chosen coordinates for the Weyl group, apart 
from explicitly h dependent terms, the c-number Hamiltonian, H{p,q), de- 
fined above, has the same functional form as the g-number Hamiltonian, 
'H{p,q). In other group coordinates that would not be the case, generally 
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speaking. Thus to associate a particular expression for a classical Hamilto- 
nian to the proper quantum Hamiltonian operator, one needs to use Cartesian 
coordinates. How can we call our choice of coordinates "Cartesian"? This 
association follows from the Fubini-Study metric induced on phase space by 
the coherent states, namely, by the fact that 

2[\\d\p,q)f-\{p,q\d\p,q)\']^dp' + dq' 

in the indicated choice of group coordinates. 

Path integrals for systems without constraints 

We now take up the question of path integrals and for pedagogical purposes 
we start with the simpler and more familiar situation in which there are no 
constraints. The abstract operator solution to Schrodinger's equation for a 
time-dependent Hamiltonian may be written as 

where T denotes time ordering. For sufficiently smooth time dependence, the 
evolution operator 

U{T) = Te^*/o 

may be represented as the limit of a large number of small time steps, namely, 
as 

U{T) = lim e""^^ • • • g-^^^^g-isWi ^ 

W— »oo 

where e = T/N, T > 0, and Hk = H{k€). This formula may be put to good 
use in at least two different ways. 

Phase space path integral — case A 

First, to form the propagator between (formal) sharp position states 
where Q\q) = q\q), for all q & R, let us insert repeated resolutions of unity 
as customary to yield 

{q"\U{T)\q')= lirn^ I ■■■ J H^o (Wil e"^^^" H^.d^^ , 
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where q" = qN+i and q' = go- As a next step we can insert resolutions of 
unity over the conjugate momentum states to give 

{q"\U{T)\q')= Xirn^j ■■■ j n^^^ (gn+i bn+1/2) (Pn+1/2 1 e-'^^"k„) 

To emphasize that we are using two different resolutions of unity, which 
requires diagonalizing both operators Q and P, and which can only be done at 
different times, i.e., sequentially, we have used the notation |g„) and \pn+i/2)- 
If we introduce the fact that 

g«Pn+l/2 9n+l 

(?n+l|Pn+l/2) = - 



27r 

as well as expand each exponential to first order in e, we are led to the 
familiar expression for the sharp q to sharp q propagator for the phase space 
path integral given by 

M j e'l^'^-^^^'^^^'^^VpVq 

-Jim y" ■ j n^^oe'^"+^/^^'"+^"'"^ [l-te{Pn+i/2\nn\qn)/{Pn+i/2\qn)] 
xn^^oO^Pn+i/2/(27r) U^=idqn . 



It is in this familiar way that meaning can be given to the formal phase 
space path integral through a close association with the abstract operator 
formulation. Of course, implicit in the expression for the overlap {qn+i\pn+i/2) 
is the assumption of Cartesian coordinates. 

While this expression is mathematically correct for a wide class of Hamil- 
tonians, it is nevertheless important to point out that it is "unnatural" from 
a physical point of view since it asserts that the phase space "paths" involved 
repeatedly oscillate between sharp q (and thereby absolutely no knowledge 
of p) and sharp p (and thereby absolutely no knowledge of q). 

Phase space path integral — case B 

We can derive another expression for the meaning of the "same" phase space 
path integral in the following way. Rather than alternately use sharp p and 
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sharp q states, let us repeatedly use just coherent states and their associated 
resolution of unity. As a consequence, the same initial expression 



= hm 6-^'"^^ • • • e-^^^2g-ieHi 



leads to 



(^",c^'\U{T)\p',q') 

= lim / • • • / n^^o (Pn+i, qn+i I e"*'^" |pn, qn) n^=i dpndqn/ {27r) , 



where p" — Pn+i, q" — qw+i and p' — Po, q' — qo- If we use the fact that 

{Pn+l,qn+l\Pn,qn) 

= exp{ii(p„+i + Pn){qn+1 - qn) - i[(Pn+l - Pnf + (Wl " ^n)^] } , 

as well as expand the exponential to first order in e, as before, we are led 
to an alternative, coherent state representation, for the formal phase space 
path integral given by 



X[l -ie{pn+l,qn+l\T-i\Pn,qn)/{Pn+l,qn+l\Pn,qn)] dpn dqn/ {2%) . 

One again, this expression is based on the implicit use of Cartesian coordi- 



Unlike case A above, this version of the phase space path integral is both 
mathematically correct for a large class of Hamiltonians as well as being 
physically "natural". It is natural because the meaning of the variables p 
and q is that of mean values rather than sharp values, and it is perfectly 
legitimate to specify the mean values of both p and q at equal times - and 
do so for all time. The meaning of these variables as mean values stems from 
the fact that {p, q\P\p,q) = p and {p,q\Q \p, q) = q. 

Remcirk: Although case A and case B led to quite different results start- 
ing from the same formal expression, it is noteworthy that they both made 





p{i{Pn+l+Pn)(«n+l-?n)/2-[(p„+l-p„)2+(9„+l-?n)2]/4} 



nates. 
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use of a first-order expansion of the exponential in the parameter e. In par- 
ticular, in both cases we made use of the approximation 

e-i^n^ ~ 1 - ieHk . 

When it comes to deal with constraints, it will become clear that this ap- 
proximation for the constraints is insufficient. 

Classical theory of constraints — a sketch 

In order to account for constraints, it is only necessary to augment the usual 
classical action functional by the addition of the constraints along with La- 
grange multipliers. The result is an action functional given generically by 
the expression 

I = Jipjq^ - Hip,q) - X"cf>^ip,q)]dt . 

Here, 1 < j < J and 1 < a < A, where J and A denote the numbers of 
canonical degrees of freedom {pj,q-^) and constraints q"), respectively, 

while A"(t) denotes the several Lagrange multipliers. Variation of Pj, , and 
A° lead to the basic equations, namely, 

qi = dH{p, q)/dpj + X"d(t)aip, q)/dpj , 
Pj = -dH{p,q)/dq^ - X''d(j)a{p,q)/dq^ , 

The subset of phase space on which the constraints holds is called the con- 
straint hypersurface. The equations of motion may also be written in terms 
of Poisson brackets. In particular, since the constraints must hold for all 
time, it is necessary that 

<i>a{p, g) = = {0a(p, q), H{p, q)} + A^{0a(p, q)Ap{p, q)} 

holds on the constraint hypersurface. This latter equation divides constraints 
into two principal classes. 

Suppose first that the Poisson brackets among the constraints vanish on 
the constraint hypersurface. In that case the second term is already zero for 
any choice of the Lagrange multipliers; the first term therefore also needs 
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to vanish on the constraint hypcrsurface (or otherwise it determines a new 
constraint that must be included). These conditions may be stated as 



Constraints that fulfill such equations are called first class constraints. A 
further division is made as follows: If the coefficients c J are constants, the 

(Jtu 

constraints are called closed first class constraints; if instead the coefficients 
c^l^ are general functions of phase space, then the constraints are called 
open first class constraints. Moreover, to solve the equations of motion it is 
generally necessary that some specific choice of the Lagrange multipliers be 
made; this is called a choice of gauge. Yang-Mills theories have closed first 
class constraints, while gravity is an open first class system. 

Next, let us suppose that the Poisson brackets of the constraints do not 
vanish on the constraint hypersurface. For simplicity, let us even assume the 
case where the Poisson brackets of the constraints {0a(p, q), <Pf3{p, q)} form an 
invertible matrix. In that case, the Lagrange multipliers are fully determined 
and are given by 



Constraints that have such properties are called second class constraints. 

Of course, there also exist mixed situations in which some of the con- 
straints are first class while the rest are second class. 

Constraint quantization — reduction 
before quantization 

In this section we outline the well known procedures of Faddeev [4] and 
Senjanovic [5] for dealing with first and second class constraint situations, 
respectively. We proceed formally as is customary in such cases. Consider 
the formal phase space path integral 



{(paip, q), <j)/3{p, q)} = cj (j)^{p, q) , 
{0a (P, H{p, q)] = hj^ (I)f3{p, q) . 



= -[{<t>a{p.q)Ap{p.q)]r' {K{p.q).H{p,q)} . 



M 





VpVqVX 
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(modulo a redefinition oi Ai), which enforces the classical constraints exactly. 
The resultant integral may well diverge (e.g., if (pi = pi and H{p, q) is inde- 
pendent of q^). Gauge fixing is used to overcome possible divergences, and 
the Faddeev-Popov determinant is introduced to maintain formal covariance 
under canonical coordinate transformations. The path integral expression 
now reads 



where x°(P)?) = 0) fo^ determines the gauge choice. This expression 

is expected to be equal to 



where B is an index that runs over the remaining, "physical" degrees of 
freedom, p* and q*. The formulation given above formally applies to the case 
of first class constraints. 

In a case of purely second class constraints, the final result is taken to be 



The foregoing expressions are plausible, formal phase space path integrals, 
but - and, in the author's opinion, this is an important qualification - these 
path integral expressions have lost any direct connection with an underlying 
abstract operator approach. While they surely can be used to calculate 
results, and on many occasions the results may well be correct, there simply 
is no firm foundation tied to an abstract operator approach to ensure that 
the results will be universally valid. 

To rectify that situation we first need to remind ourselves what is the 
accepted abstract operator formulation of quantization when constraints are 
present. 






'/^6{(l){p,q)}VpVq, 



which again is formally equivalent to an expression of the sort 
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Abstract operator quantization with 
constraints — quantization before reduction 

The general abstract operator quantization procedure for systems with con- 
straints is due to Dirac [6] . In this approach one quantizes first and reduces 
second. This is the preferred order since one then has the chance to employ 
Cartesian coordinates in the quantization, which, as described earlier, is the 
proper set of coordinates to promote to canonical operators. (Reduction first 
may give rise to a constraint hypersurface that does not admit Cartesian 
coordinates.) Thus, we suppose that we have obtained suitable canonical op- 
erators and Pj, 1 < j < J, and also chosen an acceptable factor ordering, 
if necessary, such that the Hamiltonian 7Y(P, Q) and the several constraint 
operators $a(P, Q) are self adjoint operators. Reduction consists in seeking 
a Hilbert space, S^phyg, called the physical Hilbert space, which is a subspace 
of the original Hilbert space i^, i.e., 9)phys C S). The elements of !^phys are 
those Hilbert space vectors for which 

for all a, 1 < a < A. Clearly, such vectors form a linear space. However, 
there are two special issues that must be considered. First, it follows from 
this criterion that 

[<i>a{P,Q),^p{P,Q)]%hys = 0, 

but this condition may have "^phys = as its only nontrivial solution. This 
situation arises for second class constraint systems. To deal with that, Dirac 
restricts his procedure to suitable first class systems; second class systems 
are dealt with in a completely different manner. Second, it may happen 
that the only nontrivial solutions are formal eigenvectors in the sense that 
{"^phys, "^phys) = oo. If this is the case, then some procedure must be intro- 
duced to deal with the fact that no true vectors exist that belong to S^phys- 
This procedure is not quite as straightforward as one might imagine. 

In the next section we outline a relatively new procedure [7] to deal with 
quantum constraints that is able to handle second class constraints as easily 
and with the same formalism as first class constraints, as well as having a 
well defined procedure to deal with those cases that have formal eigenvectors 
that are not in Hilbert space. 
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Projection operator method for 
quantum constraints 

Ideally, if ^a^phys — for all a, it should follow that 

holds as well. This relation works sometimes but not always. Therefore, 
let us relax this latter condition and replace it as follows. Assume that the 
operator ^a^a ^^^^ adjoint and has a spectral representation given by 

POO 

E^^l= / AcilE(A) 

expressed in terms of the associated spectral family of projection operators 
{E(A) : < A < oo}. We introduce the projection operator 

Jo 

which projects onto the spectral interval from to Here, 5{h) denotes 

a small parameter to be chosen appropriately; it is not a Dirac delta function! 
Finally, the physical Hilbert space is given by 

S:)phys = Eio . 

A few examples will help explain how the projection operator method works. 

First, let — Jk, k — 1, 2, 3, be the generators of the rotation group. 
We want to project onto those states for which Jk'^phys — for all k. We do 
so by considering 

E = M{Jl + Jl + J| < . 

Since is just the Casimir operator for the rotation group, with eigen- 

values given by j(j + j = 0, |, 1, . . ., it follows that j = is the only 
subspace allowed by the projection operator. (Clearly, a small range of other 
values for SiKf' works just as well, but we shall not dwell on that aspect.) 

Second, let $i = P and $2 = Q- The equations P'^phys = and 
Q^phys = imply that [Q,P]'ifphys = ih^phys = 0, i.e., "i/phys = 0. This 
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is the classic example of a second class system for which the original Dirac 
procedure does not work. However, let us choose 

E = E(P2 + <h) 

which acts to project onto vectors for which (Q + iP) ^phys — 0- If Q ^-iid P 
are irreducible, then the only solution is a projection onto the ground state 
of an harmonic oscillator with unit angular frequency. The essential point is 
the projection in this case is onto a one dimensional subspace. 

It is noteworthy that the first example consists of an operator with a 
discrete spectrum that contains zero (first class system), while the second 
example involves an operator with a discrete spectrum that does not include 
zero (second class system). 

Third, let $i = P be the only constraint. This operator has its zero in 
the continuous spectrum, and thus all nontrivial solutions to the equation 
P^phys — obey {^physi ^phys) — OO- In the projection operator language, 
the operator 

E = E(p2 < 5^) 

vanishes as 5 — > 0, so care must be taken to extract the "germ" of this limit. 
(An h dependence is not important in this case.) To extract the desired 
"subspace" where "P = 0" , it is most convenient to adopt a representation 
space. For that purpose let us choose a coherent state basis. In particular, 
let us consider the quotient 

(p",g"|E(P^ < 5')\p\q') /(0|E(P^ < 

= I' ^-ik-P'r/2+ikW'-q'}-ik-pr/2 ^ jT' ^-k^ 

As 5 ^ 0, the numerator and the denominator each vanish; however, the 
quotient will not vanish. Indeed, as 5 0, this quotient becomes 

which characterizes a one dimensional physical Hilbert space, which is a per- 
fectly acceptable result in this case. Since this expression no longer depends 
on q" or q', it is clear that we have reached the space where "P = 0". Ob- 
serve that the physical Hilbert space in this case is, strictly speaking, not a 
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subspace of the original Hilbcrt space Sj. Nevertheless, from a representation 
point of view, it is important to observe that the physical Hilbert space of 
interest can be obtained by a suitable limit taken from within the original 
Hilbert space Sj. 

Dynamics 

There are two important cases when dynamics is considered. The first case 
assumes that the Hamiltonian is an observable. An observable operator O is 
one which commutes with the projection operator; specifically, that 

[O, E] = . 

Therefore, if the Hamiltonian is an observable, it follows that 

[H, E] = . 

In that case, we clearly have the operator identity that 

e-^«^E = Ee-^(''^^^^E. 

This equation asserts that when Ti. is an observable and commutes with 
the projection operator, it is sufficient to impose the projection operator 
at just one time - here chosen as the initial time - and then the temporal 
evolution remains thereafter within the physical subspace, and, moreover, the 
temporal evolution is generated by that component of the Hamiltonian that 
lies within the physical subspace. The Hamiltonian is an observable for first 
class systems and for those second class systems for which the Hamiltonian 
vanishes. 

The second and more general situation is when the Hamiltonian is not an 
observable, namely, in cases for which 

[n, E] ^ . 

We would still like to ensure that the temporal evolution lies wholly within 
the physical subspace, and it is clear that one initial application of the pro- 
jection operator will not be sufficient. Just as we use the classical Lagrange 
multipliers to force the time evolving classical system back to the constraint 
hypersurface when we need to, we can use the projection operator to force 
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the time evolving quantum system back to the physical subspace when we 
need to. In symbols, this argument suggests that we consider 

lim e-''^E • • • e-*'^Ee-*'^E , 

N—^00 

where, as before, e = T/N, and T > is fixed. As shown by Chernoff [8], 
this limit is exactly 

jEg-i(EHE)Tjg ^ 

as desired. In the second class case, there are special examples where the 
temporal evolution is not unitary, e.g., liH, — P and IE is a projection onto 
the positive half line, Q > 0. However, if the original Hamiltonian is bounded 
below, which is more common in physical situations, then there is always a 
unitary version of the desired temporal evolution in the physical subspace. 

It is clear that a first class system can also be treated with repeated 
alternate projections and short time propagations, so the procedure outlined 
for second class systems works equally well for all systems. 

Integral representation for projection operator 

In special cases, such as first class systems that correspond to compact 
groups, it is straightforward to find integral representations that yield an 
appropriate projection operator. However, it it noteworthy that there exists 
a universal integral representation that yields the desired projection operator 
for any set of constraint operators [9] . We have in mind the operator identity 
given by 

which involves a time ordered functional integral over c-number Lagrange 
multipliers, where i?(A) is a suitable (weak) measure. This result holds for 
any r > (note that the left side is independent of r). The measure R{X) 
depends on r, S{hy, and the number of constraints, but it is totally inde- 
pendent of the choice of the set of constraint operators Indeed, this 
expression applies even if the constraint operators all vanish, in which case 
we learn that 
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Such an integral representation for the projection operator can be ex- 
phcitly used in forming a path integral representation for a system with 
constraints, and since the measure is the same for all systems, it may be 
used to provide a common formulation for any constrained system. Since we 
use the explicit measure in the following section, we will not describe it here. 

We now turn our attention to providing a phase space path integral for- 
mulation of temporal evolution in the presence of general constraints that 
maintains a close association with the abstract operator formulation that we 
have presented. 



Coherent state path integrals with constraints 

We wish to find an interpretation of the formal phase space path integral 
that yields the desired expression 

for temporal propagation in the physical Hilbert space. 

In the following equation chain, the weak measure -R(A) is made explicit 
as we choose a formula that achieves our goal, namely: 



■X''^aip,q)]dt 



VpVqVR{\) = lim 

AT— >oo 



N 

n 

n=0 



X lim 

M^oo 



M 

n 

m=l '- 



{Pn+m/M, Qn+m/M\Pn+{m-l)/M , Qn+{m-l)/M) 



+Sm,M{Pn+m/M, Qn+m/M\ {—ieH) \Pn+{m-l)/M, Qn+(m-l)/M) 
~\~(Pn+m/Mi (ln+m/M \ [—i 



N 



M 



X f JJ dpn+m/Mdqn+m/M j da{-fn) 

n=l ^ ■m—^ / 
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N 



n=l 

AT TV 



n=0 r-1 

lim q"\e~''^'E ■ ■ ■ e-^'^Ee-^^^E q') 



n=0 n=l 

(p",g"|Ee-^(^«^)^E|p',g') . 



Here, as usual, p",q" = Pn+i,(1n+i as well as p',q' — po,qo- The constant 
c = —AniM/e part way through the equation chain is a normalization chosen 
to ensure the form of the equation which follows the one in which c appears. 

It is important to observe that, unlike the Hamiltonian 7i, it was necessary 
to expand the expression involving the constraints to second order in the 
small parameter e. In addition, it was necessary to introduce an additional 
refinement (M) of each small time step (e) in order to construct a projection 
operator E to go along with each of the large number (N) of small time step 
evolutions for the Hamiltonian. 



Summary 

With this final expression we have achieved our goal of providing a path inte- 
gral formulation for canonical systems with general constraints that is closely 
associated with the abstract operator formulation. It is noteworthy that this 
formulation offers a path integral approach to the quantization of systems 
with first and second class constraints that does TVOT involve: gauge fixing, 
Faddeev-Popov determinants, Gribov ambiguities, moduli space, auxiliary 
variables, ghosts, indefinite metrics, Dirac brackets, etc. 
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